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ABSTRACT 


One view of the solar-wind turbulence is that the observed highly anisotropic 
fluctuations at spatial scales near the proton inertial length d p may be considered 
as Kinetic Alfven waves (KAWs). In the present paper, we show how phase¬ 
mixing of large-scale parallel propagating Alfven waves is an efficient mechanism 
for the production of KAWs at wavelengths close to d p and at large propaga¬ 
tion angle with respect to the magnetic held. Magnetohydro dynamic (MHD), 
Hall-Magnetohydrodynamic (HMHD), and hybrid Vlasov-Maxwell (HVM) simu¬ 
lations modeling the propagation of Alfven waves in inhomogeneous plasmas are 
performed. In linear regime, the role of dispersive effects is singled out by com¬ 
paring MHD and HMHD results. Fluctuations produced by phase-mixing are 
identified as KAWs through a comparison of polarization of magnetic fluctua¬ 
tions and wave group velocity with analytical linear predictions. In the nonlinear 
regime, comparison of HMHD and HVM simulations allows to point out the role 
of kinetic effects in shaping the proton distribution function. We observe gen¬ 
eration of temperature anisotropy with respect to the local magnetic held and 
production of field-aligned beams. The regions where the proton distribution 
function highly departs from thermal equilibrium are located inside the shear 
layers, where the KAWs are excited, this suggesting that the distortions of the 
proton distribution are driven by a resonant interaction of protons with KAW 
fluctuations. Our results are relevant in configurations where magnetic held inho¬ 
mogeneities are present, as, for example, in the solar corona where the presence 
of Alfven waves has been ascertained. 


Subject headings: 
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Introduction 


Turbulence in plasmas is a complex phenomenon which is characterized by different 
regimes in different ranges of spatial and temporal scales. Turbulence in the solar wind 


has been extensively studied, both b y det ailed ana. 


a theoretical point of view; see 


yses of in-situ measurements and from 


Bruno fe Carbonel (120051 1 for a review. Such studies often 


adopt complementary views that the turbulence may be described either as collection of 
wave that interact nonlincarly, so-called wave-turbulence, or else as a collection of broad 


band, essentially zero frequency eddies or flux tubes that form_a hierarc 


structures. These approaches have been extensively reviewed (Barnes 


ly of coherent. 


1979; 


Matthaeus et al. 


2015) and we do not attempt a critical comparison in the present work. Instead, we adopt 


mainly a wave taxonomy of the fluctuations, based on linear theory in order to address 
a specific se t of questions. As motivation, we not e that a variety of observations in the 


solar wind (IBale et al. 


2005 


Sahraoui et ah 


2009) have suggested that fluctuations near 


the end of the magnetohydrodynamics inertial cascade range, and approaching the kinetic 
plasma range, may consist primarily of Kinetic Alfven waves (KAWs). Here we address in 
particular the nature of fluctuations produced due to nonlinear interactions near the proton 
inertial length d p and investigate in some detail the basis for identifying them as KAWs. 
We show how phase-mixing of large-scale parallel propagating Alfven waves is an efficient 
mechanism for the production of KAWs at wavelengths close to d p and at large propagation 
angle with respect to the magnetic field. To support the interpretation as KAWs, we 
perform and analyze Magnetohydrodynamic (MHD), Hall-Magnetohydrodynamic (HMHD), 
and hybrid Vlasov-Maxwell (HVM) simulations that model the propagation of Alfven waves 
and their fully nonlinear interaction with a nonuniform plasma background. We will be 
able to characterize fluctuations produced by this “phase-mixing-” like interaction as highly 
oblique KAWs. 
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At frequencies much lower than the proton gyrofrequency, solar wind turbulence is 
dominated by Alfvenic fluctuations, which are characterized by highly-correlated velocity 


and magnetic fields and by low-level (with respect t o t 


and magnetic field intensity variations flBelcher fe Davis 


re background values) density 


197111 . Moreover, in recent years 


several indications have been found of the presenc e of Alfven waves also in the solar corona 


( To mczyk et ah 


20071 : 


Tomczyk & McInto sh 


2009 ) from where the solar wind emanates. 


In the solar wind, Alfvenic fluctuations extend over a wide range of spatial scales, from 
large scales within the MHD range, down to lengths comparable or smaller than the proton 
inertial length d p , where also finite Larmor radius effects become relevant. Moreover, in this 
scenario it is important to take into account the role of spectral anisotropy. Theoretical 
studies, indeed, have shown that in a MHD turbulence the energy c ascade p referentially 


takes p lace perpenc 


Carbone fe Veltril (1990); 


icula rl y to the background m agnetic field B 0 (e.g., 


Shebalin et ah 


Oughton et al 


4198,1k 


(1994)). Then, it is expected that at smaller 


scales the fluctuation energy tends to concentrate in nearly perpendicular wavevectors. 
This idea is supported by observational data showing that the distribution of wavevectors 


of magnetic fluctuations h as a s i 


magnetic field flMatthaeus et ah 


gnificant population quasi-perpendicular to the ambient 


1986 . 


1990). All the above effects, such as turbulence 


in a plasma with several characteristic lengthscales, Alfvenic correlations and spectral 
anisotropy, can interact with large scale structures such as shears and magnetic equilibria, 
with a subsequent generation of wave-like activity. The full picture need to be addressed 
with plasma simplified models, as well as the Vlasov-Maxwcll treatment. 


Waves belonging to the Alfven branch, with wavelengths comparable with the proton 
inertial length d p and wavevectors nearly perpendicular to the mean magnetic field B 0 
are often indicated as “kinetic Alfven waves”. During the last decades, KAWs have 
received considerable attention and have been studied in detail due to their possible role 
in a wave description of the turbulent cascade. Since the MHD cascade favors nearly 

































5 


perpendicular wavevectors, the expectation within a wave perspective would be that 
KAWs are naturally p resent at sc ales of the order of d p . An extensive analysis of KAW 


physics is found in 


Hollwed (119991) (see also references th erein for a more 


on the subject). Many solar wind observational analyses (IBale et al. 


20091 : 


Podesta & Tenbarge 


theoretical works (Howes et al. 


2012 ; 


as numerical simulations uGarv fe Nishimura 


Salem et al 


2008a: 


2012 


Chen et al. 


2004: 


2013 


2005 


complete view 


Sahraoui et ah 


Kivani et al. 


2013), 


et al. 

2009] 

Sahraoui et al. 

2012 

) as wcl 

Howes et al. 

2008b; 

TenBarge & Howes 


2012) have suggested that KAWs can play an important role in the dissipation of 
turbulent energy. Due to a nonvanishing electric held parallel component associated with 


KAWs, these waves have a 

so been considered in the Droblem o 

particle acceleration 

(Voitenko & Goossens 

2004 

Decamp & Malara 

2006 

). Recently, 

Vasconez et al. 

(2014) 


have studied collisionless Landau damping and wave-particle resonant interactions in 
KAWs. 

There are also simplified problems, less complex than fully developed turbulence, in 
which one finds the formation of small scales in the direction perpendicular to an applied 
magnetic held B 0 . It is well known that this effect appears in the context of MHD when 
imposed parallel propagating waves interact with an inhomogeneous ba ckgrou nd consisting 


either of pressure balanced structures or velocity shears ( Ghosh et al. 


1 9 98). In 2D 


equilibria, where the Alfven velocity varies in directions perpendicular to the magnetic held, 


two mechanisms have been investigated in detail: (1) phase-mixing (Heyvaerts & Priest 


19831), in which differences in group velocity at different locations progressively bend 


wavefronts; and (2) resonant absorption which concentrates the wave energy in a narrow 
layer where the wave frequency locally matches a characteristic frequency (Alfven or 
cusp). These processes h ave been studied both by investigating normal modes of the 


inhomogeneous structure fl Kappr aff fc Tataronis 


1985 


Davit 


1987 


Hollweg 


1987 


Califano et ah 


977; 


Mok & Einaudi 


1990 


1985 


Steinolfson 


19921 ) and by considering the 
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evolution of an initial disturbance (Lee & Roberts 


1986; 

Malara et al. 

1992, 

1996) 


of density stratification and magnetic line divergen ce (lRuderman e t al 


nonlinear coupling with comp ressive moc 


localized pulses flTsiklauri fc Nakariakov 


es (Nakariakov et al. 


2002 


Tsiklauri et al.. 


997 


998|), as well as 


1 9981 ), and evolution of 


2003) have been considered. 


The propagation of MHD wave s in inhomogeneou s magneti c fields containing null p oints 


has also been studied in detail (jL an di et al. 


( 20051) ; see also 


McLaughlin et al. 


feoiol) for a 


review). Phase-m ixing in 3D inhomogene ous equilibria has also been co nsidered in the smal 


wavelength limit (Similon & Sudan 


Malara et al 


19891) using a WKB approxi mation (Petka 


deta 


20001) . also within the problem of coronal heating flMalara et al 


2003 


1998 


2005 


20071 ). Par ticle acceleration in ph ase-mixing of Alfven waves in a dispersive regi me has bee n 


studied by 


Tsiklauri et al. 


2011) and in 3D (Tsiklauri 


( 2005 ) using particle-in-cell simulations, both in 2D flTsiklauri 
20121) configurations. Finally, instabilities generating KAWs in 


plasma with transverse density modulations have been considered by 


Wn fc Chenl (120131 ). 


Similar ideas involving dissipative mechanisms related to interaction of Alfven waves or 


KAWs and phase-mixing 


lave been examined in 


sheet (ILvsak & Song 


he cont ext o f the mag netosp heric plasma 


2011) and in coronal loops flOfman & Aschwanden 2002). 


The above considerations suggest that phase-mixing of Alfven waves might represent 
a mechanism to produce KAWs, when the wavelength of waves becomes enough small to 
be comparable with the proton inertial length d p . This effect could work, for instance, in 
the solar corona where the background magnetic field is clearly inhomogeneous and where 
the presence of Alfven waves has been ascertained. This represents an indication about the 
nature of small-scale fluctuations that could be present in the coronal plasma. The aim of 
the present paper is to investigate directly the generation of KAW fluctuations associated 
with Alfven wave phase-mixing. This study will be performed numerically by using both 
a HMHD code and a Vlasov-hybrid code. The former include proton skin depth effects, 
while the latter allows identification of kinetic effects such as wave-particle resonances 
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and temperature anisotropy. In Section 2 the model will be presented along with a linear 
analysis of wave properties; the results derived by the HMHD in the linear regime are 
described in Section 3, while the nonlinear regime and the results of the Vlasov-hybrid code 
will be described in Section 4; a discussion and a summary of results are given in Section 5. 


2. Hall-MHD model 


Consider a fluid plasma composed of protons and electrons. For phenomena at 
sufficiently large scales the electron mass can be neglected m e m p and the quasi¬ 
neutrality condition holds: n e ~ n p , with n e and n p the number densities of electrons 
and protons, respectively. Moreover, we assume that the temperatures of protons and 
electrons are equal T e = T p = T. Finally, the displacement current term is neglected in the 
Ampere equation. In these conditions the plasma dynamics can be described by the HMHD 
equations: 


<9v 

dt 


| + V.(pv)=0 

(v • V)v = ——V(pT) + - [(V x B) x B] 
2 P P 


SB 

~dt 

dT 

dt 


= V x 


vxB- (VxB)xB 

P 


+ (v- V)T+(7-l)T(V- v) = 0 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


Equations (JH)-(j4j) contain only dimensionless quantities. In particular, mass density p 
(which is only due to protons) is normalized to a typical density p, temperature T to a 
typical value T, the pressure p = pT due to both protons and electrons is normalized to 
the value p = ‘IkbPT/ m pi with kb the Boltzmann constant and m p the proton mass. The 
spatial coordinates are normalized to a typical length L , magnetic held B is normalized to a 
typical magnetic held B, fluid velocity v to the typical Alfven speed ca = B /(47T/7) 1 / 2 , time 
t to the Alfven time = L/ca■ Finally, the plasma beta [3 = p/(B 2 /8n) is a typical value 




for the kinetic to magnetic pressure ratio; 7 = 5/3 is the adiabatic index; e = d p /L -C 1 
is the Hall parameter measuring the relative amplitude of the Hall term with respect to 
the v x B term in the Ohm’s law, d p = Ca/^cp — CA'm p c/(qB) being the proton inertial 
length. Equations (fT]) - (Hj) reduce to the ordinary single fluid compressible MHD equations 
in the limit e —> 0. The equations are ideal and adiabatic due to the absence of viscosity 
in Eq. (2), resistivity in Eq. (3) and heat sources and heat conduction in Eq. (4). For the 
simulations carried out below these effects are present in some form (not necessarily in the 
form of a fluid model); however at present they are omitted for clarity. 


2.1. Wave properties 

We consider Erst a uniform equilibrium state characterized by homogeneous 
dimensionless density po, temperature T 0 , magnetic held B 0 and vanishing velocity v 0 = 0. 
Waves propagating in the above equilibrium can be studied as usual: a small amplitude 
perturbation is superposed on the equilibrium; equations ([I])-®) are linearized with respect 
to the perturbation amplitude and perturbation fields are Fourier transformed both in space 
and time. Imposing nonvanishing perturbations leads to the following dispersion relation: 

w 6 + CV + C 2 u 2 + C 3 = 0 (5) 

with 

c ao^n ( k f \ + k l ) 

Ci = -q (24 + 4) - 4 (4 + 4)-Wj!- >- ro 

(fcjj + 2 k^kj + k\ 

^p0 

c 3 = -c 4 A 0 c 2 s 0 kf\ (kl + kl) (8) 

In equations, (JS])-® uj is the wave frequency normalized to ; kn and k _l are the 
wavevector components parallel and perpendicular to B 0 , respectively, both normalized to 


C 2 = 


r 2 k 2 
c A0 fi \\ 


C A0 + 2c so) { k \\ + k l) + 


c 4 r 2 k 2 
c A0 c s0 ,i \\ 
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Z -1 ; cao = -B 0 /(47rpo) 1 ^ 2 and c s o = (/ S7T0/2) 1 / 2 are the Alfven speed and the sound speed 
associated with the equilibrium, respectively; co p0 = ( qBB 0 /(m p c))(L/cA ) is the normalized 
proton gyrofrequency. The squared sound speed to Alfven speed ratio will be indicated by 
A) = c 2 s0 /c\ 0 . 


The expres sions 


Vasconez et al.l (2 014 ). i 


m-m are equivalent to those found in a two-fluid model by 


solved (Vasconez et al 


the electron inertia is neglected. Equation (J5]) has been analytically 


20141) using the Vieta’s substitution method (IBirkhoff fc Mac Lane 


19771 ) for the calculation of the complex roots of a third-degree algebraic equation. For given 
values of kn and k± the three solutions found for co 2 are real and positive and correspond to 
the Alfven, fast magnetosonic (FM) and slow magnetosonic (SM) branches, respectively. I 11 
particular, KAWs correspond to the Alfven branch for k\\ -C k±_ ~ co p q/cao- 


From the linear analysis the expressions for the amplitudes of fluctuations can be 
calculated. We use a reference frame where the x-axis is along Bo while the wavevector k is 
in the xy plane. The perturbation components for a wave with wavevector k = k\\e x + k±e y 
can be expressed in the following form: 

t ’ 1 * = V (X- - c y?) - 4 ^ ) sm{k " x + kxv - ut + 4i) (9) 


'AO 


A 11 ^pO 


V\ y = a- 


CO 


2 02 l, 2 

Vo^ii 


CO 


co 2 (co 2 — c 2 0 k 2 ) 


02 k 2 


sin(A;||a; + k±y — cot + 0) 


'AO 


viz = -a- 


c A0 k \\ 


CO 


cos(fc||a: + k±y — cot + 0) 


B lx = aB { 


Bi y - 


h j j k±uo p0 


co 


cok 2 

-aBM (l - 


02 1,2 
-ifvii 


sin (A; 11 re + kj_y — cot + (ft) 


-A0 r 


CO 


cok 2 




02 1,2 

\rv\\ 


sin(fc||X + k\_y — cot + 0) 


'A0 f 


Biz = aB 0 cos(fc||£ + k±y — cot + (ft) 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 
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where a is the perturbation amplitude, </> E [0, 27r] is the phase and lu = cj(k) is derived 
from the dispersion relation (131) . 


In what follows it will be useful to have an expression for the component of the 
wave group velocity perpendicular to B 0 , v g j_ = duj/dk j_. This has been calculated in 
the following way: we indicate by F(u,k\\,k±) the LHS of equation (J3D and by oj(k\\,k±) 
a solution of equation (j3J). Substituting co(k\\, k±) in the place of u into equation (J5]) we 
obtain an identity 

F(u(k\\, k±), k\\, k±) = 0 (15) 


which holds for any value of k\\ and k±. The derivative of equation (1131) with respect to kj 

is 

dF dF duo dF 

+ =0 


from which we obtain 


dk _|_ du dk i_ dkj_ 

du dF/dk± 

V 9 ±- ~ 


dkj_ dF/ du; 

The derivatives of F appearing in the RHS of equation (fT6l) can be calculated using 
equations (|5])-(f8j); their expressions are: 


(16) 


dF 


dk 


= -2 




k± (4 0 + c^ 0 ) + 


k ± kfA 0 


UJ, 


pO 


UJ A + 


k±k\\ (2c 2 s0 c A0 + 4 0 ) + 2 


k±kt + k\kl ) c\ 0 c% 


UJ, 


pO 


(U - 2k±knC sO c A0 (17) 


dF 

dcu 


= 6u 5 - 4 


( C A0 + C s0 ) + 4 c A0 + 


/V / ..2 


UJ', 


pO 


cu 3 + 


k 2 k\ (2 c 2 s0 c 2 A0 + 4 ) + 


/.2l,4 4 2 

"'ll"' L UO c sO 


CU 


pO 


UJ 


(18) 


where k 2 = k 2 ^ + Ar[. Since phase-mixing increases A;_l leaving fc|| constant (as shown in the 
following section), in Fig. Q]the perpendicular group velocity v g ± of the Alfven, FM and SM 
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modes is plotted as a function of k± in a range k± > k\\, for k\\ = 1 and for (3o = 1.25 and 
/3q = 2.08. The choice of these two values of /3 0 will be justified in the next section. From 


Fig. [T]we see that v g ± of the three modes increases with increasing k±, becoming nearly 
constant at large propagation angles. For k± > ku — 1, v g ± is positive for Alfven and FM 
waves while is negative for SM waves. Moreover, v g ± of FM waves is much larger than that 
of the Alfven branch; this is consistent with a quasi-isotropic propagation velocity of FM 
waves. The behavior of v g ± for negative k±_ can be inferred from Fig. |Tj taking into account 
that v g ±(k\\, — k±) = — v g ±(k\\, k±). These features will be used to identify the nature of 
fluctuations generated by phase-mixing. 


3. MHD and Hall-MHD simulations of phase-mixing 


We consider a 2.5D configuration where all the physical quantities depend only on two 
spatial variables (x and y ), but vector quantities can have three nonvanishing components. 
The above quantities are defined in a spatial domain D = {(x,y)} = [ 0 , 27 t ] x [0, 27 t ] (in 
dimensionless units), where periodic boundary conditions are imposed both in the x and y 
directions. We consider a nonuniform equilibrium structure in D, where physical quantities 
vary only along the ^-direction, or are uniform. Quantities relative to the equilibrium are 
indicated by the upper index ”(0)”. The equilibrium magnetic field is 


B (0) = B (0) ( y) 


(19) 


where 



( 20 ) 
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and e x is the unit vector along x. The dimensionless parameters in the expression 
the following values: b m = 1.5, h = 0.2, r = 10 and 

(b m - l)r 


a = 


2(2 h) r 


1 + 


1 

2 h 


r-i 2 


~ 2.62 x 10 


-4 


have 


( 21 ) 


The function B^\y) is symmetrical with respect to the central point y = ir, where it 
reaches its maximum value b m . The small term containing a in the expression (120|) has 
been added in order to have a vanishing first derivative of B^(y) at the two boundaries 
y — 0 and y = 27 T. Both B^\y) and dB^/dy are periodic functions in the interval [0, 27 t]. 
However, higher order derivatives of B^°\y) are not exactly periodic; as a consequence, the 
Fourier spectrum of B^\y) has a tail at high wavenumbers. To avoid this drawback, the 
expression ( 1201 ) has been corrected by filtering out harmonics with wavenumbers larger than 
70 in its spectrum. The filter does not sensibly alter the profile B^\y). The equilibrium 
magnetic held is almost constant both in the central half and on the two lateral parts of the 
domain, while two sharp shear layers (current sheets) are located in between these uniform 
regions. The equilibrium temperature, which is equal both for protons and electrons, has 
been chosen as uniform, while the fluid velocity is identically zero: 


T (°) = x ; v (o) = o 


( 22 ) 


The equilibrium mass density p^(y) is determined by total pressure equilibrium: 

^\v)t^ + ^3M =P m (23) 

where P^° = 1.748 and /3 = 2 . The Alfven and sound velocities associated with equilibrium 
structure are given by c^\y) = B^°\y)/[p < '°\y)} 1 ^ 2 and c® = (/TyT^/2) 1 / 2 , respectively. 
The local plasma (3 is ^°\y) = [cf^/c'^y)} 2 . The profiles of c^\y) and /3^(y) are shown 
in Fig. [2] The Alfven velocity is larger in the lateral parts of D than in the center; the 
inhomogeneity of is responsible for phase-mixing of Alfven waves. We note that (3^ < 1 
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in the central region, while f3^ becomes larger than 1 when approaching the boundaries 
y — 0 and y = 27T. In particular, /3^ ] = 1.25 in the middle of the shear layers while 
/ 3 = 2.08 in the lateral homogeneous regions; these two values have been used to calculate 
the profiles of Fig.sQ]and[5] 

At the initial time an Alfvenic perturbation has been superposed on the above 
equilibrium. Quantities relative to the perturbation are indicated by the upper index ”(1)”. 
The initial magnetic held and velocity perturbation are given by 


B (1 ^(a;, y, t — 0) = a cos(x)e z 


v (1) (a;, y,t = 0) = -a[p (0) (y)] 1/2 cos(x)e z (24) 


so that v*9) = -(c®/B®)BW. The quantity a gives the amplitude of the initial 
perturbation. Initial density and temperature fluctuations are vanishing: p^(x, y, t = 0) = 
0, T<U(x,y,t = 0) = 0. 

The MHD case is described by equations (jTT) -(141) with e = 0. In this case and in the 
small-amplitude limit a 1 the equations may be linearized and above initial conditions 
evolve in time according to the equations: 


B (l \x,y,t) = a cos x-c^(y)t e z , v {1) (x,y,t) = -a[p (0) (y)] 1/2 cos x~c { ^\y)t 


( 0 ), 


(25) 


indicating that the initial perturbation propagates along B (0 ) at the local Alfven speed. 
Thus, transverse variations of c® generates in the perturbation increasingly smaller scales 
which are localized within the shear layers. This phenomenon represents phase-mixing of 


an Alfven wave (IHevvaerts fc Priest 


1983). 


The fully nonlinear equations (fl jl -fpf jl have been solved numerically with the above- 
specified boundary and initial conditions. The HMHD numerical code employs a 2D Fourier 
pseudospectral method to calculate spatial derivatives and time integration is performed via 
a second-order Runge-Kutta scheme. Aliasing errors in the evaluation of nonlinear terms 
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are partially removed by a 2/3 truncation in the spectral space. Equations have been solved 
using the scheme described in 


Ghosh et ahl (119931 ). Hyper-viscosity and hyper-resistivity 
terms (fourth-order derivatives) have been added in equations ([2]) and (J2D in order to obtain 


numerical stability with dissipation concentrated only at the smallest spatial scales. 


Different runs have been performed; Table 1 summarizes the values of parameters used 
in the various cases, with n x and n y the number of gridpoints in the x and y direction. 
RUN 1 corresponds to the MHD case (e = 0) with a low-amplitude perturbation (a = 0.01). 
Though the analytical solution for an infinitesimal amplitude is known (equations ()25l) ). we 
performed this run to single out differences of a purely MHD case with respect to a HMHD 
case (e / 0). Fluctuating fields are defined as Sg = g — ( g ) x , where g represents a physical 
quantity from which the spatial average along x has been subtracted in order to eliminate 
the contribution from the equilibrium structure. The effect of phase-mixing is visible in 
the time evolution of 8v z and 5B Z (not shown): the wave propagates from left to right but 
with a velocity which is larger in the central part than in the lateral parts of the domain. 
As a consequence, within the two shear layers the wave profile is stretched with small-scale 
gradients generated in the transverse ( y ) direction. The other components of <5v and 5B, 
as well as the density and temperature fluctuations, have a much lower amplitude which is 
of the order of 10~ 4 , or smaller. Moreover, these fluctuations have a wavelength in the x 
direction which is half of the wavelength of 5v z and 8B Z . These two features clearly indicate 
that in the MHD case Sv x , 8v y , 8B X , 8B y , 8p and 8T are generated by small nonlinear 
effects that are quadratic in the Alfven wave amplitude. 

In RUN 2 we considered again a small-amplitude initial wave (a = 0.01), but now 
dispersive effects have been switched on by setting the coefficient of the Hall term e = 0.125. 
This corresponds to set the proton inertial length d p = 0.125 (in normalized units). Thus, 
it is expected that dispersive effects are no longer negligible as soon as the perturbation 
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wavevector k has increased enough (in consequence of phase-mixing) to become comparable 
with dp 1 . This condition is reached inside the shear layers. The time td at which k ~ d p 1 


can be estimated in the following way: during phase-mixing t 
component k y increases in time according to the equation (e.g. 


re transverse wavevector 


Petkaki et ah 


(1199811 1 


k x t (26) 

d y ) s 

where ( dcAo/dy) s is an estimation of the Alfven velocity gradient in the shear layers, while 
k x is constant and k y0 — k y (t — 0). In our case k x — 1, k y0 = 0 and ( dcAo/dy) s — —1 (in 
the shear layer at x ~ 4.5). At the time t = t d we have k y k x , then, from the condition 
ky(td)d p ~ 1 we get the estimate 


ky(t ) ~ kyQ 


td - 


1 


cLL 


f dcAo\ 


~ 8 


(27) 


" *V Ay J, 

In our configuration it is expected that at times t >t d phase-mixing could generate KAWs, 


i.e., perturbations belonging to the Alfven branch with a quasi-perpendicular wavevector 
(k y k x ). We now explore this possibility in some detail. 


In Fig. [3] the fluctuating fields Sp, ST, 5v and <5B are represented at the time t = 13.1, 
along with the current density component j z which has been plotted in order to localize 
the equilibrium shear layers. The time t — 13.1 is larger than t d , thus we expect to observe 
effects due to finite ion inertial length. We observe that Sv x , 5B X , Sp and ST are now of 
the same order as Sv z and SB Z and the parallel wavelength is 2tt for all these fields. This 
indicates that in the HMHD case Sv x , SB X , Sp and ST fluctuations are not due to nonlinear 
effects but they are part of the same perturbation as Sv z and 5B~, namely, a KAW. A 
clearer identification of this perturbation as a KAW will be given in the following. 


In Section 2.1 we found that the perpendicular group velocity of KAWs is nonvanishing, 
though it is much smaller than the background Alfven velocity. Then, in the present 
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configuration, KAWs generated inside the shear layers, while propagating along B (0 \ slowly 
drift in the y direction. The perpendicular group velocity v g ± of KAWs has the same sign 
as k y = k _l (see Fig. [TJ, the latter being negative (positive) in the shear layer located 
at y ~ 1.8 (y ~ 4.5). Then, KAWs produced in both shear layers would move outside 
toward the lateral higher- (3^ uniform regions. Indeed, examining the time behavior of the 
perturbation, oblique wavefronts progressively occupying the two lateral regions are found. 
This behavior can be seen in Fig. [4] where the 5B X component is plotted at four different 
times. We used the following procedure to measure the propagation of these structures in 
the direction perpendicular to B^ 0) . 


Once the system undergoes phase-mixing, and once wavepackets are generated (at 


about t ~ 3), we identify all the local maxima and minima of t 


fields give similar results). In practice, following 


Donato et ah 


re density field p (different 


(119981 ). all the critical points 


x* where Vp = 0 have been found. In these points we computed the square Hessian matrix 
of p, identifying the strongest maxima and minima. We followed the trajectories of these 
points in time (which are almost parallel each other), that go from the shear layer out 
to the border. Taking the average position between the maximum and the minimum, we 
calculated its velocity, obtaining u y =~ ±5.4 x 1CD 2 . We can also estimate their wavevector 
components k y = ±2n/\ y ~ ±8.11 and k x = 1 . Using these values in the expression (TTB1) 
we obtain for the perpendicular group velocity of the Alfven branch (see Fig.s|Tjand[ 2 ]) the 
value v g± = ±7.1 x 1(T 2 for 0(°) = 2.08 (in the lateral homogeneous region), in reasonable 
agreement with the value u y estimated in the simulation. We conclude that the observed 
structures propagate with the group velocity of KAWs. We note also that the propagation 
angle is 6 = arctan \k y /k x \ ~ 83°, close to 7 t/2 as required for KAWs. 


From Fig. |T]we note that v g ^_ is opposite to kj_ for waves belonging to the SM branch. 


This implies that SM waves possibly produced in the shear layers would laterally drift 
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opposite to KAWs. Thus, the waves observed in the lateral homogeneous region cannot 
belong to the SM mode. Finally, Fig. [T] indicates that for large propagation angles the 
group velocity of FM waves is much larger than that of KAWs, being ~ 1.8 for the above 
values of k x and k y . This value is much larger than the lateral propagation velocity u y of 
the observed perturbation. Consequently, the fluctuation produced by phase-mixing cannot 
belong to the FM mode. 

Another feature we took into account to identify the observed perturbation is the 
polarization. Since the wavevector of these perturbations is nearly parallel to the y 
direction, the condition V • B = 0 implies that the dominant components are Bi x and B\ z . 
We considered these components in the polarization analysis. Equations (jT2|) and Q 
indicate that the magnetic perturbation is elliptically polarized. In particular, using the 
dispersion relation ([5]), we calculated the quantity 

^(k, A)) = 1 T~l2 ( 28 ) 

c ao k \\ 

appearing in equations (TT2l) and (fT4lh In Fig. [5] we report a as a function of k± for the three 
modes, for k\\ — 1 and for two values of (5 o, corresponding to the lateral uniform region 
(/3o = 2.08) and to the middle of the shear layer (f3 0 = 1.25). From Fig. [5] we see that a 
is negative for Alfven and FM waves, while is positive for SM waves. Then, for positive 
k\\ and k± and given values of t and x, as one increases the y coordinate, the perturbation 
magnetic field rotates clockwise (counterclockwise) in the zx plane for Alfven and FM 
waves (SM waves). This characteristic behavior of the eigenmodes has been compared with 
the simulation results. 

In Fig. [6] we plot a hodogram in the SB X - 5B Z plane, parameterized by the coordinate y 
which varies in the range \n, 2n\. This sample is for a fixed value of the x-coordinate x = i r, 
and at the time t — 20. At that time, the perturbation generated inside the shear layer 
almost fill the lateral uniform region. In the hodogram the blue diamonds indicate the shear 
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layer (4 < y < 5) while red triangles indicate the lateral uniform region (5 < y < 2n). The 
blue asterisk indicates the central point y = tt and the red square indicates the boundary 
y = 27r. It is seen that in most of the shear layer and in the lateral homogeneous region 
the perturbed magnetic field clockwise turns with increasing y. Since in this region k y 
is positive, the observed polarization in the simulation is in accordance with that of the 
Alfven branch. 

In conclusion, based both on the group velocity and on the polarization analysis, 
we deduce that the fluctuations generated inside the shear layers are KAWs. From the 
hodogram we also notice that the perturbation amplitude in the shear layer is smaller than 
in the lateral uniform region. Then, the process of KAW generation and their subsequent 
lateral propagation tends to move the initial Alfven wave energy away from the shear 
regions. A process of local fluctuating energy depletion is eventually found also in the pure 
MHD case, but it is simply due to dissipation localized at the shear layers. 

In the hodogram of Fig. [6] a small perturbation can be seen at the boundary 
between the shear layer and the inner homogeneous region, in which the magnetic field 
counterclockwise turns. This can be interpreted as a SM perturbation propagating from 
the shear layer toward the central region, i.e., opposite to k y: as predicted by the linear 
theory (see Fig. [[]). This SM perturbation can be seen also in Fig. [4] as a fluctuation 
drifting toward the center y = n of the domain. However, the energy associated with this 
perturbation is much lower than that of the fluctuation that we have identified as a KAW. 


4. Large-amplitude HMHD and kinetic simulations 

Next, we consider the phase-mixing of an initial large-amplitude (a = 0.25) Alfven 
wave. We show results from both a HMHD (RUN 3) and a kinetic (RUN 4) simulation. 
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The results are qualitatively similar as in the small-amplitude case, as it can be see in 
Fig. [7] where 5B x (x,y ) is illustrated in the x, y-plane for the four runs at time t = 13.1. 

In particular, also in the large-amplitude case phase-mixing of the initial wave generates 
small scale variations perpendicular to B (0 ^ in the shear layers, mainly in the form of 
KAWs. The identification of these waves employs the same method as in the low-amplitude 
case, namely, considering both the perpendicular group velocity of perturbations and their 
polarization. The main difference between the low and high-amplitude case is that in the 
latter a small-amplitude precursor of the main KAW perturbation is observed to fill the 
lateral homogeneous region before the arrival of the main perturbation (Fig. [71). 


ation (RUN 4) has been performed using a HVM numerical code 
. The HVM algorithm integrates numerically the Vlasov equation for 
the proton distribution function in multi-dimensional phase space. In the present work, we 
restrict our analysis to the 2D-3V (two dimensions in physical space and three dimensions 
in velocity space) phase space configuration. The electrons are considered as a fluid and 
a generalized Ohm equation is employed for computing the electric field, which retains 
the Hall term. In the present work electron inertia effects are neglected. Quasi neutrality 
is assumed and the displacement current is neglected in the Ampere equation, therefore 
assuming low frequency dynamics. Finally, an isothermal equation of state for a scalar 
electron pressure is employed to close the HVM system. The equations solved by the HVM 
code are the following: 


The kinetic 

simu 

(Valentini et ah. 

2007 


| + u.V/ + i(E + uxB).| = 0 
at e <9u 


E = —v x B + — [jxB — -VF f 

n ° 


<9B 

~dt 


-V x E 


V x B = j 


(29) 

(30) 

(31) 


where f(x,y,u x ,u y ,u z ,t) is the proton distribution function in phase space and E (x,y,t) is 
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the electric field. The proton density n and the ion bulk velocity v are obtained as velocity 
moments of /. The scalar electron pressure P e is derived from an isothermal equation 
of state assuming that the electron temperature is equal to the initial (uniform) proton 
temperature. In equations (l29lh(l3Tl) all quantities are dimensionless, as specified for the 
HMHD equations (JTJ)-((4]) ; moreover, the velocity u is normalized to ca, the density n to 
p/m p , the electric field E to caB/c, and the current density j to cB/(AtiL). A detailed 
description of the numerical method employed to solve equations (T29|) - d3T]) can be found 


m 


Valent ini et ah, 


(2007J). The 2D-3V phase space has been discretized (see table 1) with 
n x x riy = 256 x 1024 grid points in the spatial domain and 51 3 grid points in the velocity 


domain. 


The initial and boundary conditions (in the physical space) used in RUN 4 are the 
same as those used in the previous large-amplitude HMHD run (RUN 3). In the 3D velocity 
domain, the distribution function / is set equal to zero for |u| > u max , where u max = 5 v t h, P 
and i>th, P is the proton thermal speed. Moreover, the initial unperturbed proton distribution 
function is a Maxwellian with a uniform temperature T^ 0) = 0.5. From Fig. [7] one can 
observe that the fluctuations 5B X obtained from the HVM run (RUN 4) at time t = 13.1 
are similar to those seen in the HMHD run (RUN 3): KAWs develop in the shear layers 
and slowly drift toward the lateral uniform region. This similarity is probably related to 
the fact that both models (HMHD and HVM) employ the same form for the Ohm’s law. 
In Fig. [8l we compare the power spectra of In \SB\ 2 from RUN 3 (left panel) and RUN 4 
(right panel), at the time t = 13.1. The contour plots in this figure display similar features, 
i.e. marked anisotropy along the k y direction due to the phase-mixing process, but we can 
notice that more energy has transferred into small scales in RUN 4 relative to RUN 3. This 
suggests the presence of enhanced small-scale activity when kinetic effects are retained in 
the description of the plasma dynamics. 
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The ion microscopic dynamics which is described by the HVM model introduces 
evidently new effects with respect to the HMHD model, which are described in the following. 
First, we notice that the amplitude of 8v x at t = 13.1 is much lower in RUN 4 than in RUN 
3. This difference can be better appreciated in Fig. |9l where profiles of Sv x at as a function 
of y at x = 7r and at two different times are plotted. It can be seen that the amplitude of 
8v x , which is initially vanishing, increases in time in the HMHD run, while it remains at a 
much lower level in the HVM run. Such a different behavior is presumably due to kinetic 
damping effects in RUN 4 that act on velocity fluctuations parallel to the background 
magnetic field; such damping mechanisms are absent in HMHD. 


The HVM models allows us to follow how the proton distribution function is distorted, 
due to resonant interaction of protons with the KAW fluctuations, at different space 
positions and times with respect to the initial Maxwellian (in contrast, the HMHD model 
assumes local thermodynamic equilibrium). To characterize the departure of the computed 


f from a Maxwellian distribution, we define an L 2 -norm difference flGreco et al. 


2012 


Valentini et al. 


2014): 


e(x,y,t) 



[/(x, u, t ) — M (x, u, t )] 2 d 3 u 


(32) 


where M is a Maxwellian with the same density, bulk velocity and isotropic temperature as 
/. e is a positive definite quantity and may be viewed as a “distance” or separation between 
the computed / and an equivalent Maxwellian. In Fig. [TO] the maximum 


£max(t) = max e(x,y,t) (33) 

(x,y)eD 

is plotted as a function of time. We see that £ max increases in time, eventually saturating 
at a value ~ 0.035; this indicate a progressive departure from a Maxwellian distribution. 
In Fig. [TO] the time t c i is indicated as a vertical red-dashed line, corresponding to the time 
necessary for phase-mixing to produce transverse wavevectors comparable with d~ l . At 
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t ~ td the growth of £ max becomes slower, indicating that the largest departure from a 
Maxwellian is almost reached when fluctuations at scales of the order of the proton inertial 
length are formed. This is an indication that the KAW fluctuations are responsible for 
the modifications in the proton distribution function. Another quantity describing the 


departure of / from a Maxwellian is the temperature anisotropy parameter ( Perrone et al. 
2013 ) 


R(x,y,t) = 1 - 


T p ±(x,y,t) 


(34) 


T p \\(x,y,t) 

where T p ± and T p \\ are the proton temperature perpendicular and parallel to the local 
magnetic field B, respectively; R < 0 (R > 0) corresponds to T p j_ larger (smaller) than 
T p ||. From results of RUN 4 we found that e and \R\ are clearly correlated, the correlation 
coefficient being ~ 0.74. This indicates that the departures of / from a Maxwellian 
are essentially due to the generation of temperature anisotropy. In Fig. Qj]the spatial 
distributions of R (left) and e (middle) are plotted at the time t = 13.1, together with 
the parallel electric field (right); comparing with Fig. [7] we see that the largest departures 
from a Maxwellian are spatially correlated with KAWs, i.e., with fluctuations at transverse 
scales comparable with d p . Positive and negative variations of the anisotropy parameter R 
with similar amplitudes follow one another along the KAW profile. Then, a prevalence of 
parallel or perpendicular proton kinetic energy can equally takes place, due to interactions 
with KAWs, according to the wave phase. A 3D representation (surface plot) of the proton 
velocity distribution at time t — 13.1 is shown in the top row of Fig. [T21 for two different 
spatial locations: (x,y) = (5.4,1.5) (left), where R is minimum at the given time, i.e., 

T p l > T p \\, and (x,y) = (6.2,4.7) (right), where £ is maximum and T p ± < T p \\. As it is 
clear from these two plots, both distributions depart from a spherical shape, typical of a 
Maxwellian distribution, and display a structuring in the form of rings perpendicular to 
the local magnetic field, indicating groups of particles in resonance with large-amplitude 
fluctuations. In the bottom row of the same figure, the contour plots (together with the level 
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lines) of the proton velocity distribution in the u x -u y plane (for u z = 0) are presented at the 
same time and spatial locations as for the top row plots. In particular, the bottom-right 
panel, corresponding to the location where the ratio T p \\/T pl _ is particularly large, shows 
the presence of a well-defined beam of ions moving in the direction parallel to the local 
magnetic field. Such a beam is spatially localized where the KAW fluctuations have been 
generated by phase-mixing, and therefore it is reasonable to postulate that the beam is 
produced by a wave-particle interaction between protons and the fluctuations that we have 
characterized as KAWs. 


Such interaction could be related to the fluctuating parallel electric field 5E\\ associated 
with the KAW. In fact, the Ohm’s law (1301) allows for the presence of an electric field 
component E\\ = E • B /B parallel to the magnetic field. This component is due to the 
electron pressure gradient term in equation (1501) . In Fig. fill E\\ is plotted at time t = 13.1. 
It can be seen that a fluctuation 5E\\ ~ 10~ 2 forms mainly in the shear layers, where 
KAWs are localized. The presence of parallel electric field fluctuations is another feature 
characterizing KAWs. In the ideal HMHD case Eu has no effects on the dynamics (the 
term containing VP e /n being canceled out when calculating VxE), unless one includes 
also a resistivity. However, in the HVM model it can have an influence on the evolution 
of the proton distribution function. In order to investigate this possibility, we considered 
the potential energy variation (per particle) 5U associated with 5E\\, in comparison with 
the proton thermal energy E t h■ The ratio between these two quantities, expressed in our 
normalized units, can be estimated as: 


5U 


8 SE iiAi 


E, 


th 3 e/3 T p 


(35) 


where Am is the parallel wavelength and T p is the proton temperature. Using the values 


e = 0.125, j3 = 2, SEu ~ 10 , Am = 27t, and T p ~ 0.5, from equation 


we find 


SU / E t h ~ I- Then, the potential energy variation associated with the parallel electric field 
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associated with the KAWs is comparable to the proton thermal energy. This implies that 
5E\\ is able to sensibly modify the initial proton distribution function. Moreover, the velocity 
of the proton beam is Ub eam ~ 1.8. We measured the phase velocity of the perturbation along 
x, finding w x ~ 1.6 which is comparable with Ubeam■ These considerations strongly suggest 
that the observed field-aligned particle beam is generated by the resonant interaction of 
protons with the parallel electric field 5E\\ associated with the KAW fluctuations. 


5. Discussion and conclusions 


In this paper we have shown that fluctuations having the character of oblique Kinetic 
Alfven Waves are readily generated using a simple configuration consisting of an in-plane 
two dimensional sheared magnetic field, and an out-of-plane perturbation that locally 
propagates as an Alfven wave. We have described the emergence of the KAWs as occurring 
due to phase-mixing, or refraction of the perturbation wave vectors towards angles highly 
oblique with r espect to t he sheared magnetic field direction. Although the refraction is 


seen in MHD flGhosh et al. 


19981) . the basic physical picture of the emergence of KAWs 


requires at least a model as complete as HMHD, where it is seen at both low and high 
initial perturbation amplitudes when the thickness of perturbation across the magnetic field 
becomes comparable with the proton skin depth. For the case of a Hybrid Vlasov solution, 
at this stage additional features are observed during this process, such as non-Maxwellian 
proton distributions, temperature anisotropy, and the formation a parallel beam in the 
proton velocity distribution. Evidence for resonant wave-particle interaction between the 
beam and the KAW perturbation is also identified. 


The simplicity of the initial setup adopted in these numerical experiments made 
it possible to characterize features of KAWs in several ways. These included, for 
example, examination of wave polarizations and phase speeds based on linear theory. The 
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identification was facilitated by the choice of the magnetic shear as an ideal MHD, ffMffD 
and ffVM equilibrium, as well as the choice of perturbation as an Alfven mode. However, 
it is equally important to recognize that much of the physical picture described in this 
idealized context is expected to carry over to more complex configuration, including even 
a full quasi-incompressible turbulence cascade. To elaborate on this point briefly we note 
that incompressible couplings, such as those that dominate the present models, correspond 
to quadratic terms, for example in the MHD equations (Eqs (2) and (3)) written for 
constant density. These nonlinear couplings (we are speaking here of the equations prior to 
separation into equilibrium and perturbation) have the familiar property that the Fourier 
amplitude with wave vector k may interact directly with two other wave vectors p and q 
provided that the triadic condition k = p + q is satisfied. Suppose we identify the sheared 
magnetic field in Eq. (19) as = (A?f°)( 2 /),0,0) having Fourier modes at wavevectors 
p —> (0,p y ,0). Now let the perturbation consist of wave vectors q and require that it be of 
an incompressible type but otherwise relax the restrictions of the special case given in Eq. 
(24). It is immediately clear that even in a large amplitude situation more complex than 
what we treated above, the triadic nonlinear couplings will drive excitations to wavevectors 
k = (0,^,0) + (q x ,Qy,Qz) that will acquire increasing values of k y . That is the fluctuations 


will spectrally transfer towards wavevectors that are oblique to 


field direction. This is closely related to the standard argument flShebalin et ah 


he sheared m agneti c 


1 983 ) for 


perpendicular spectral transfer, and is very similar to the nonlinear phase-mixing associated 
with velocity shears and pressure bala nce s tructures seen in MHD turbulence simulations 


with analogous initial setups ( Ghosh et ah, 


19981 ). I 11 this way we can see that the driver of 


the refraction towards highly oblique wave vectors is considerably more general than in the 
special case we considered here. It is reasonable to suppose then, that the Hall and kinetic 
effects that we identified as emerging when the transverse scales decrease to the ion inertial 
length will also be observed in these more complex circumstances. 
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We can conclude then that the simple mechanism we described here in the context 
of HMHD and Hybrid Vlasov models provides a pathway to understand generation of 
fluctuations with the character of Kinetic Alfven Waves. Driven by magnetic shear, 
the excitations appear at ever smaller scales across the magnetic shear layer until the 
Hall/kinetic effects appear near the ion inertial length scale. Analysis of polarization and 
propagation speed identifies these fluctuations as being of the KAW-type. Particularly 
interesting is the emergence of a field-aligned beam in the self-consistent proton distribution 
function, apparently admitting signatures of wave-particle resonances. We conjecture that 
the mechanism described here operates also in a quasi-incompressible cascade scenario, so 
that the basic reasoning given here may also account for emergence of KAW-like features 


(Bale et ah 

2005; 

Sahraoui et al. 

2012) 
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Fig. 1.— The perpendicular group velocity v g ± of the Alfven (black lines), FM (red lines) 
and SM (blue lines) inodes is plotted as a function of k±, for k\\ = 1 and for j3o = 2.08 (full 
lines) and /3q = 1.25 (dashed lines). 
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Fig. 2.— The Alfven velocity (black solid line) and /3^ (red dashed line) associated 
with the equilibrium structure are plotted as functions of the y coordinate. 
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Fig. 3.— Sp (a), ST (b), j z (c), 8v x (d), 5v y (e), Sv z (f), 8B X (g), 5B y (h), SB Z (i), plotted 
as functions of x and y, at the time t = 13.1 for the low-amplitude HMHD run (RUN 2). 
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Fig. 4.— The 5B X component is plotted in the xy-plane at four different times for the 
low-amplitude HMHD run (RUN 2). t = 2 (a), t — 10 (b), t = 15 (c), t — 20 (d). 
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Fig. 5.— The quantity a is plotted as a function of k±, for = 1 and for /3 0 = 2.08 (full 
lines) and /?o = 1-25 (dashed lines). Black lines correspond to the Alfven mode, red lines to 
the FM mode, and blue lines to the SM mode. In the latter case the curves corresponding 
to the two values of f3 0 are superposed. 
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Fig. 6.— The 5B Z and 5B X components are plotted on the two axes for y varying in the 
range n < y < 2i r, at x = x 0 = vr and t = 20. The central point y = tt (blue asterisk), 
the boundary y = 2tt (red square), the shear layer region (blue diamonds), and the lateral 
homogeneous regions (red triangles) are indicated. 
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Fig. 7.— The component 8B X is plotted as function of x and y, at the time t = 13.1, for (a) 
RUN 1 - MHD with small a; (b) RUN 2 - HMHD with small a; (c) RUN 3 - HMHD, large 
a; and (d) RUN 4 - HVM large a. 
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Fig. 8.— Power spectrum of In \ 5B\ 2 from RUN 3 (left panel) and RUN 4 (right panel) at 
t = 13.1. 
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Fig. 9.— 5v x profiles are plotted as functions of y, for x = tt, at times t — 1 (left panel) and 
t = 13.1 (right panel), for RUN 3 (black full line) and RUN 4 (red dashed line). 
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Fig. 10.— The quantity £ m axi quantifying the departures of the HVM distribution function 
from a Maxwellian, is plotted as function of the time. The vertical dashed line indicates the 
time t = td, described in the text. 
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Fig. 11.— The quantities R (left panel), e (middle panel), and the electric held parallel 
component E\\ (right panel) are plotted in the au/-plane, at time t = 13.1. 
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Fig. 12.— Top row: surface plot of the proton velocity distribution at the spatial location 
(x,y) = (6.2, 4.7) (left), where e is maximum and T p ± < T p \\, and (x,y) = (5.4,1.5) (right) 
where R is minimum, i.e., T p ± > T p \\. The magenta tubes in the two plots indicate the 
direction of the local magnetic field. Bottom row: shaded contours (together with level 
lines) of the proton velocity distribution, in the u x -u y plane (at u z — 0), in the same spatial 
locations as in the top row. 
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Table 1: Simulations setup. 


RUN 

Type 

Spatial Resolution (n x x n y ) 

Amplitude (a) 

Hall parameter (e) 

1 

MHD 

256 x 256 

0.01 

0 

2 

HMHD 

256 x 256 

0.01 

0.125 

3 

HMHD 

256 x 1024 

0.25 

0.125 

4 

HVM 

256 x 1024 

0.25 

0.125 








